Continuous fourier transform method and apparatus

ABSTRACT

An input analog signal to be frequency analyzed is separated into N number of simultaneous analog signal components each identical to the original but delayed relative to the original by a successively larger time delay. The separated and delayed analog components are combined together in a suitable number of adders and attenuators in accordance with at least one component product of the continuous Fourier transform and analog signal matrices to separate the analog input signal into at least one of its continuous analog frequency components of bandwidth 1/N times the bandwidth of the original input signal. Given the separated frequency components, the original analog input signal can be reconstituted by combining the separate analog frequency components in accordance with the component products of the continuous Fourier transform and analog frequency component matrices. The continuous Fourier transformation is useful for spectrum analysis, filtering, transfer function synthesis, and communications.

United States Patent [191 Munoz Nov. 26, 1974 [75] Inventor: Robert M. Munoz, Los Altos, Calif.

[73] Assignee: The United States of America as represented by the Administration of the National Aeronautics and Space Administration, Washington, DC.

[22] Filed: Apr. 18, 1973 [21] Appl. No.: 352,382

[52] US. Cl. 235/156, 235/197, 324/77 B [51] Int. Cl. ..L G06f 15/34 [58] Field of Search 235/156, 152, 184, 181, 235/197 [56] References Cited UNITED STATES PATENTS 3,344,349 9/1967 Schroeder 324/77 3,529,142 9/1970 Robertson 235/197 X 3,610,901 10/1971 Lynch 235/156 X 3,636,333 1/1972 Klund 235/156 3,714,566 1/1973 Kang 235/156 X T910009 5/1973 Abramson 235/156 X [57] ABSTRACT An input analog signal to be frequency analyzed is separated into N number of simultaneous analog signal components each identical to the original but delayed relative to the original by a successively larger time delay. The separated and delayed analog components are combined together in a suitable number of adders and attenuators in accordance with at least one component product of the continuous Fourier transform and analog signal matrices to separate the analog input signal into at least one of its continuous analog frequency components of bandwidth l/N times the bandwidth of the original input signal. Given the separated frequency components, the original analog input signal can be reconstituted by combining the separate analog frequency components in accordance with the component products of the continuous Fourier transform and analog frequency component matrices. The continuous Fourier transformation is useful for spectrum analysis, filtering, transfer function synthesis, and communications.

22 Claims, 10 Drawing Figures /5 V a e" f fil l T N 4 12 r an) DEI AY N t /3 i A DElT AY I at) 1 1 DELAY I PATENTELHUYZB'B" Y 3,851.1 e22- sum MIF 4 TRANSMITTING RECEIVING END END 2/ CHANNEL 1 CHANNEL 2 INVERSE ON IN OUS INPUT L F CONTINUOUS OUTPUT m MECHANISM TRANSFORM f (T) MECHANISM CHANNEL N V 24 INVERSE A CONT'NUOUS CONTINUOUS FOUR'ER FOURIER O TRANSFORMER T R C (T) 22 .coN i Uo US FOURIER-TRANSFORM M T-non 1 AND APPARATUS l The invention described herein was made an employee .of the United States Government and maybe manufactured and-used by or for the Government for governmental purposes withoutthe'payment' of any royalties thereon or therefor. r

BACKGROUND OF THE INVENTION The present invention relates general .to' method and apparatus for performing Fourier transformations .ing discrete frequency components in the frequency domain, thereby obtaining a discrete spectral representation of the input signal to be analyzed. It is also'known from the prior art that once an analp'g signal (has been. separated into its v corresponding discrete Fourier components in the frequency domain that a'discrete representation'of the originalanalog'signal can-be reconstituted by performing an inverse Fourier transformation on the discrete spectral Fourier components in a 'digital'computer to derive a discrete repr'e sentation'of the 'original analog signal in thetime domain. The discrete'fast Fourier transform'met'hod and apparatus is disclosed in an article What is the Fast Fourier Transform?, appearing in the I.E.E.E. Trans- I matrix (D matriIQ'and the Fourier components matrix V .analogsignal are combined accordance with the component productlof the discrete Fourier transform (C,,,(t) matrix);to signal.

reconstitute the original analog input In another feature of the present invention, an'analog input signal having a certain bandwidth is Fourier transformed into N number of Fourier frequency components each component having l/N X'the bandwidth of the original-input signal. Such separate Fourier fre- "quency components are transmitted simultaneously over N number of different channels to a receiver. The

actions on "Audio and Electroacoustic's,"volume I AU-15,No. 2, of June 1967, pa e 45-55. a

On'e of the problems withthe discrete Fourier trans form isthat. it requires the use of 'a digital computer with associated analog t'oj-digit'al convertersfand the 3 analysis performed by the computer comprises'many compl'eir multiplications, additions, ,and memory cycle -M SUMMARY OF THE PRESENT INVENTION The principal object of the present invention is the provision of an improved Fourier transform method and apparatus and more particularly of a simplified I continuous Fourier transfonn method which does not require a complex digital computer for implementa- I tion. I

' In one feature of the present invention, ananalog sig-- a. na] to be analyzedis separated into N number of analog components each delayed relative to the original'signal or the reference one'of said signal components by a successivly larger time delay, such reference and de- I layed components are added together in accordance with at ieastone of the component products of the fast Fourier-matrix (D matrix and the delayed input signal FIG. 2,

' quency,

received Fourier components are recombined in accordance with the inverse Fourier transformto reconstitute the 'originalinput signal, whereby a relatively wide bandwidth signal istransmittedoverN number of relatively narrowband individual channels. v

Other features and advantages'of the present tion will become apparent upon a perusal of the following specification taken in connection .with the accompanying drawingswher'ein:

I BRIEF Esc'R PrIoN. OF THE DRAWINGS FIG. 1 is a signal flow graph depicting the continuous Fourier transformation of a complex (two-wire) analog input signal into four separate complex Fourier frequency components,

FIG. 2 is a signal flow graph depicting the continuous Fourier transformation of a real. (one-wire) input signal to its respective Fourier frequency components,

FIG. 3 is a physical implementation of the system of FIG. is an alternative physical implementation of the input portion of the circuit of FIG. 3 and delineated by line 44, I

' FIG. Sis a plot of; the inputsignal to be Fourier analyzed and the output Fourier components as a function of time and depicting the discrete Fourier transform output indicated by-the dots superimposed upon the continuous Fouriertransform output,

FIG.' 6 is a schematic circuit diagram for a physical implementation of 'a continuoius Fourier transform cir- I similar toQthat' FIG; 3 exceptfor eight time deh i put n. p

FIG." 7 is a'composite plot of band limited input signal data, frequency bands of the separate Fourier components, and the'frequency response of the first Fourier frequency output coefficient C ,all as a function of fre- FIG. 8 is a schematic block diagram of a data transmission system incorporatingfeatures of the present invention, I v

FIG. 91is a diagram similar-to that of FIG. 8 for modifying an input signaL'and FIG.(10 is a schematic-circuit diagram of an inverse continuous- Fourier transform apparatus.

, DESCRIPTION OF THE PREFERRED EMBODIMENTS.

.- "The prior method ofdiscrete'Fourier transformwill first be described followed .by a description of an extencomponent matrix F(t to produce at least one continuous output function C, ,,(r) corresponding to one of the frequency components of the analog input signal robe analyzed. I I

In another feature of the present invention, the Fou rie'r componentsin the frequency domain of an original sion thereof under thepresent invention to a continu- 'us' Fourier transform method and apparatus. In the discrete Fourier-transform method, a band limited continuous complex function of time fl'tlis consideredto contain .j a maximum frequency component co /2 :wherein w'gis'the rate-atwhich thefunction'flt) is to be invensampled to produce a sampled data function f*(nT), where T l/w seconds. A set of N sequentially ordered magnitudes off*(nT) where n O, l, N l, is selected and the discrete Fourier transform of this set is as follows:

t A 1 N DFT (f( ))=cm= 2 e where m=O, 1,. .Nl

and c are the N complex frequencies coefficients.

The inverse operation of computing the original data set from the coefficients is given as follows:

21rymn f* m 1;) N Eq. 2

CNI(I) where n=0, l, N-l

These equations can be written in matrix form in the following way:

' and the fast Fourier transform reduces this by a factor of 2 log N/N which is appreciable for large N.

In the continuous transform method of the present invention, the discrete Fourier transform equations (3) and( 4) are utilized andthe input signal f(t) is considered as a vector quantity having real and imaginary components. Because Equations (3) and (4) are independent of time, the independent variable F can be made a function of time F (t) to produce the following continuous Fourier transform:

C(t)=1/NDF(!).

Matrix Equation (5) may be expanded and written as follows: i

m) gfufl' As in the fast Fourier transform a considerable sim- C l/ND F (3) plification in the computation is achieved if N is selected to be a binary number, i.e., a power of 2 such as and 2, 4, 8, 16, etc. This simplifies F and greatly reduces the complexity. The matrix Equation (6), when solved for F =D C the complex frequency coefficients C,,,(t), produces a q. set of equations as follows:

00) =l/Nlfo( )+f|(l) 'l'fzU) "fv-|( )l i(l)=l/Nlfl1(l)+ fi( f2(l) fiv UH C2(t) =l/Nm(i)f,(i) +f (1) g -"f- ,(r)| 9) iv-i( )=l/ lfi:( e" fv-i(r)l (1 where C l l l e e: gli j For N 4, these equations reduced to: f e e e 0( LIOU) +f1( +f2( +f3( )l q- (11) 6" D tN-HJ g undF= C (r)= l/N muwaflupnupefion Eqv(l2) f l T) 2( li) The discrete fast Fourier transform is obtained for C t 1 N t Equations (3) and (4) by recognizing and using the 3( [M emu) f2) +8730 EH14) properties of the D matrix to simplify the computations. Since all entries in the D matrix are periodic The signal flow graph for Equations 11-14 is illusroots of unity, and symmetric about the main diagonal, trated in FIG. 1. In the signal flow graph, each node a great simplification is possible where N is a power of represents a variable, and the arrows terminating at 2. The highest degree of symmetry and hence the greatthat node originate at the nodes whose variables conest simplification is permitted when N fulfills this retribute to the value of the variable at that node. The quirement. The matrix equations of (3) and (4) are contributions are additive, and the weight of each contypically solved with a digital computer programmed to tribution, if other than unity, is indicated by the conperform the fast Fourier transform. The program instant written close to the arrowhead of the transmisvolves many complex multiplications and additions. Al-

though the computations are relatively complex the sion. Equations 1 l-l4 are separated into real and imagmary components as follows:

mm Rem) re ulat s 11 Equations 15-22 could be implemented with hardware, however, considering a one dimensional input;

signal having only a real and no imaginary component,

Equations 15-22 are greatly simplified by recognizing. that:

ReC (t) ReC (t) l/N [Refl,(t) Ref t)] Eq. (25) Since ImC (t) and ImC (t) zero, no hardware implementation for Fourier frequency components is necessary or possible. Therefore, equations 23-30 reduced to only four independent equations, namely Equations I3, 25, 26 and 27, the signal flow graph for which is depicted in FIG. 2. In the continuous transform, f (t) is identical to the reference input signal f(t) and fl(t) is identified with f (t+T), f (t) is identified with fi,(t+2T) and f (t) is identified with f (t+3 T).

The signal flow graph of FIG. 2 is implemented with hardware as shown in FIG. 3 to provide a continuous Fourier transform apparatus. More particularly, the input function of time flt) (magnitude) is first divided by N in attenuator 11 to provide a first reference input signal fl,(t)/N. UN is a scaling constant equal to the reciprocal of the number of Fourier frequency components. A sample of the output of the attenuator 11 is fed to a first delay 12 to be delayed by a time T. The output of the delay 12 forms the second input signal, namely f (t)/N. A sample of the output of delay 12 is fed to a second delay 13 which delays the signal passing therethrough by the same time T such that the output of the second delay 13 fonns the third input signal, namely f (t)/N. An output of the second delay 13 is fed to a third delay 14 for delaying the fourth signal by a time Trelative to the third signal to produce the fourth input reference signals namely f (t)N The continuous input signals are combined in the inverting operational amplifiers 15 in accordance with the signal flow graph of FIG. 2 to produce the four outputs, namely, ReC0(t), ReC2(t), ReC (t) ReCaU), ImC1(t) One relatively simple apparatus for deriving the time delayed input signals for the input to the circuit of FIG. 3 is as shown in FIG. 4. More particularly, the input signal flt) to be Fourier analyzed is recorded on a magnetic recording tape 16. The magnetic recording tape 16 is drawn from a supply reel 17 onto a take-up reel 18 across a number of pick-up heads 19 each head 19 being equally spaced from the adjacent head in the di rection of movement of the recording tape 16. The outputs from each of the downstream heads 19 is successively delayed by a successively larger time compared to the output of the upstream head or reference head. The respective output signals are indicated above each of the respective heads.

Referring now to FIG. 5, there is shown the signal to be analyzed fl!) and its Fourier frequency components for N 4 and corresponding to the four outputs of the circuit of FIG. 3. In FIG. 5,flt) is shown being separated into N 4 components and sampled four times per time window where the time window No. I is defined as O to 4T. The information contained in the four components ReC (t), ReC (t), ReC (t) and 1mC (t) is no greater nor less than that contained in the four original samples because the continuous Fourier transform process is conservative, that is, it permits information neither to be created nor destroyed. Since only one unique set of Fourier frequency components per window is allowed to characterize flt), according to the sampling theorem, each Fourier frequency component must then contain only l/N X the bandwidth limit of flt or a maximum frequency content of (u /2N (u /8 Hz. This accounts for the bandwidth reduction in the discrete or sampled data case and this has a direct and immediate implication for the continuous transform mechanism.

The continuous transform method is related to the series of discrete Fourier transforms in successive time windows 1, 2, p and at the limit of this process, i.e., when p goes to infinity and T goes to zero, the bandwidth relationships are invariant for both 1 and T. Therefore, they obtain in the limit and apply equally well for the continuous transform case, i.e., the bandwidth of any Fourier component C (t) l /N bandlimit of f(t).

The continuous Fourier transform mechanism of the present invention is suited to many practical applications including filtering, spectral analysis, and transfer function modification in the frequency or time domain. In short, it can do all that the theory permits of any Fourier transform mechanism and for some applications it results in very simple and inexpensive hardware.

Referring now to FIG. '6 there is shown a Fourier transform circuit similar to that of FIG. 3 where N 8. As can be seen in FIG. 6, the hardware implementation involves merely attenuators and amplifiers as contrasted with the prior art discrete Fourier transform mechanism that requires complex multiplication, addition and memory cycles. The frequency characteristics of the continuous transform mechanism for the circuit of FIG. 6 is shown in FIG. 7 where N 8 and Imf(t) O. The information band offlt) from zero to (u /2, is separated into five equibandwidth parts by the output frequency components C through C, as indicated by frequency bands 0 through 4. This operation of producing Fourier frequency components on a continuous basis is identical to the operation of filtering with the added advantage that an exact measurement of the phase of the incoming signal with respect to w, is possible for each band. It can be seen that the C component acts as a low pass filter and C acts as a bandpass filter, etc., each band having a bandwidth of (+)m,/2N around their respective center frequencies of mw /N where m O, l (N1). All this filtering action takes place simultaneously and continuously and therefore constitutes a parallel set of M2 filters or a spec trum analyzer. If only one filter is needed, then the hardware can be simplified greatly.

At frequencies above (0 /2 the action of filtering, unlike that of ordinary filters, repeats itself in a cyclical manner as indicated by the frequency band containing the same number designation in FIG. 7. This was the quency components. The components are then transmitted over separate channels 1 through N to be reconstituted at the receiving end via an inverse continuous transform mechanism 22 into a facsimile of the original signal.

Referring now to FIG. 9 there is shown an apparatus 23 adapted for transfer function modification in the frequency domain. More particularly, the system is similar to that of FIG. 8 with the exception that variable attenuators 24 are placed in each of the respective channels between the continuous Fourier transform mechanism 21 and the inverse continuous Fourier transform mechanism 22 for modifying the frequency components in the frequency domain. The reconstituted modified output signal G(t) appears at the output across a load resistor 25. Convolutions in frequency or time can be performed by using two continuous transform systems of the type shown in FIG. 9 and multiplying the input signals or the frequency components respectively and reconstituting or transforming these products using another transform mechanism.

Once an input signal in the time domain has been separated into its continuous fast Fourier transform components Cm(t) in the frequency domain these components can be recombined to form the original signal f(t) by means of the inverse Fourier transform, the discrete version of which is shown in Eq. (4). In the continuous Fourier transform, Equation (4) is rewritten as:

F(t) D C (t) reason for band limiting the input signal flt) to avoid 25 Equation (31) is solved for fi (t) =flt) and for N 4 signal aliasing and the consequent confusion of data as follows:

13(1) fli) 1 i 1 C t!) f,(r) flt-T) l 2 e C,(1)

Eq. (32) f n) flr2T) l e 0* A f 1) f( I3 T) l e e C3( 1) that results. However, for some applications this action From Equation (32): may be desirable for use in a periodic filter as, for exflz) C (r) C,(r) C (4) C (t) Eq. (33) ample in the magnetic amplifier or flux gate magm- Taking into account real and imaginary components tometer where the output signals are contained in all even harmonics of the drive signal and it is necessary flt) Rec) ImCOm Ream Imam to reject the odd harmonics and extraneous noises in Reczfl) ImRzm Recam ImRsm (34) other bahds- Assuming the imaginary part of F t) 0 then, lmC (l) Many types of periodic filters are possible such as 1mC2( t) O and Rea) Recs) and 1 g even and odd harmonics as well as third, fifth, eleventh, ]mC3( then; thlrteenth, nineteenth harmonic types, each seflt) c 2 ReClU) R8620) (35) lected as combinatlons of frequency bands containing in the same band number or numbers. This action can, Equahoh 15 Implemented y the hardware of of course, be accomplished simultaneously for different 10 and Cohslsts y of the Positive Sum of three combinations and may be useful for human speech and pohehts' The amphfier 22 is a hon-inverting p music Signa] analysis tional amplifier with four input terminals each of a Returning now to band limited information input sig- 5O lf gaihinput Signal u) is connected lhto nals, the continuous transform mechanism can be used two 1hPut termlhals to Yleld a galh of The Output to k m i i System bl f di of the inverse continuous transformer 22 is the reconhigh bandwidth information over a plurality of N low stituted Original input Signal fi bandwidth lines or channels. Such a transmission sys- The continuous Fourier transform mechanism of e tem i Shown i FIG 8, I FIG 8, h input i l j(;) present invention is useful for filtering, signal process is continuously Fourier transformed in Fourier transg and Communications In filtering, it y be utilized former 21 in a manner as previously described with refo OWPaSS filtering, bandpass filtering, parallel filtergard to FIG. 3. The output of transform mechanism 21 ing and chorus filtering (a selected set of harmonically is N frequency components with the attendent bandrelated frequencies and combinations), specialized filwidth reduction of l/N for each of such Fourier fretering of spectral signatures, and phase lock filtering or phase discrimination such as that used in space communication, and it is also useful for complex signal filtering.

The advantage of the continuous Fourier transform apparatus of the present invention as applied to filtering is that very accurate frequency selection is obtained there being no practical limit on the minimum bandwidth. The filtering circuits are inexpensive to implement in hardware as the hardware comprises merely amplifiers and in some cases attenuators.

In signal processing, such as in a system described with regard to FIG. 9, the continuous Fourier transform apparatus of the present invention is useful for transfer function modification in the frequency domain. It may also be utilized for spectrum analysis of continuous data (quasi stationary and non stationary time series produce deterministic results). It may also be utilized for serial-to-parallel data conversion and parallel-toserial data conversion.

In communications, such as the system described with regard to FIG. 8, the continuous transform mechanism is useful for bandwidth reduction in the individual communication channels. It may also be useful for video tape recording utilizing multiple tape heads.

What is claimed is: 1. In a method for continuous analog analysis for a time varying input function flt) in the time domain into at least one of its Fourier Frequency components C,,,(t) the steps of:

continuously producing N signal components f (t)/N ranging from f (t)/N to f ,(t)/N wherein each signal component except f (t)/N is delayed an amount T from the previous signal component;

continuously combining separated ones of said signal components [f,,,(t)] f (t)/N together so as to produce at least one continuous output function C,,,(t) corresponding to one of the Fourier frequency components of the time varying input function flt) to be frequency analyzed.

2. The method of claim 1 wherein the step of combining said separated components f (t)/N together so as to produce at least one output function C,,.(t) corresponding to one of the Fourier frequency components of the time varying input signal flt) comprises the step of, combining the time delayed continuous signal components f (t)/N together according to at least one of the component products of the matrix product: DF (t) where D is the discrete Fourier transform matrix and F (t) is the delayed input signal component matrix.

3. The method of claim 2 wherein the step of combining the delayed signal components f (t)/ N according to at least one component product of the matrix product includes the step of, adding the delayed signal components f,,,(t)/N according to at least one of the component products of said matrix product.

4. The method of claim 2 wherein D is defined by the matrix:

5. The method of claim ll wherein N is a power of 2.

fi D where D is the discrete Fourier transform matrix and C(t) is the matrix for the continuous separate Fourier frequency components.

7. The method of claim 6 wherein the step of combining a plurality of the continuous time varying Fourier frequency components C,,,(t) together so as to reconstitute the original time varying function flt) includes, combining the separate Fourier frequency components C,,,(t) according to the product of the matrix Equation:

where D is the discrete Fourier transform matrix and C(t) is the matrix for the continuous separate Fourier frequency components.

8. The method of claim 7 wherein the D matrix is defined by the matrix equation:

10. In an analog apparatus for continuous frequency analysis of a time varying input function flt) in the time domain into at least one of its Fourier frequency components C,,,(t):

means for continuously generating N signal components f,,.(t)/N ranging from f (t)/N to f (t)/N wherein each signal component except fi,(t)/N is delayed an amount T from the previous signal component;

means for continuously combining said signal components f (t)/N together so as to produce at least one continuous output function C,,,(t) corresponding to one of the Fourier frequency components of the time varying input function flt) to be frequency analyzed.

11. The analog apparatus of claim 10 wherein the means for combining said separated signal components f,,,(t)/N together so as to produce at least one output function C (t) corresponding to one of the Fourier frequency components of the time varying input signal flt) includes, means for combining the time delayed continuous signal components f,,,(t)/N according to at least one of the component products of the matrix product:

DF(t) where D is the discrete Fourier transform matrix and F (t) is the separated input signal component matrix.

12. The analog apparatus of claim 11 wherein said means for combining the delayed signal components according to at least one component product of the matrix product includes, means for adding the delayed signal components according to at least one of the component products of said matrix product.

13. The analog apparatus of claim 11 where D is defined by the matrix:

14. The analog apparatus of claim wherein N is a power of 2.

15. Apparatus comprising:

first analog means for continuously producing N time varying Fourier frequency components C t) in the frequency domain of a time varying function flt) in the time domain;

said first analog means comprising: means for continuously generating N signal components f (t)/N ranging from f (t)/N tof (t)/N wherein each signal component except f (t)/N is delayed an amount T from the previous signal component; and means for continuously combining said signal components f,,,(t)/N together so as to produce N varying Fourier frequency components C,,,(t);

channel means for conveying each time varying Fourier frequency component to an inverse Fourier transformer, said transformer comprising second analog means for combining the continuously time varying Fourier components C,,,(l) together so as to reconstitute the original time varying function flt).

16. The apparatus of claim wherein said means for combining a plurality of the continuously time varying frequency components C,,,(t) together so as to reconstitute the original time varying function flt) includes, means for combining the separate frequency components C,,,(t) according to the component products of the matrix product DC(t) where D is the discrete Fourier transform matrix and C(t) is the matrix for the continuous separate Fourier 18. The apparatus of claim 17 wherein:

C00) 10) C(z) 19. In anan alog apparatus for continuous frequency analysis of a time varying input function flt) in the time domain into at least one of its Fourier frequency components C,,,(t):

means for continuously generating N signal compo nents f,,,(t)/N ranging from fl (t)/N to f (t)/N wherein each signal component except f (t)/N is delayed an amount Tfrom the previous signal component;

said generating means including a tape recorder with a recording head and (N-2) staggered playback heads, said signal component fi,(t)/N being connected to said recording head and each of the remaining signal components being produced at one i of said playback heads;

means for continuously combining said signal components f,,,(t)/N together so as to produce at least one continuous output function C (t) corresponding to one of the Fourier frequency components of the time varying input function f(t) to be frequency analyzed.

20. An analog apparatus as set forth in claim 19 wherein said means for continuously combining signal components consists of a network of operational amplifiers.

21. Apparatus comprising:

first analog means for continuously producing N time varying Fourier frequency components C,,,( t) in the frequency domain of a time varying function flt) in thetime domain; y

said first analog means comprising: means for continuously generating N signal components f,,,(l)/N ranging from fi t)/N to f- (t)/N wherein each signal exceptf (t)/N is delayed an amount T from the previous signal component; and means for continuously combining said signal components f,,,(t)/N together so as to produce N varying Fourier frequency components C,,,(t); said generating means including a tape recorder with a recording head and (N-2) staggered playback heads, said signal component fl,(-t)/N being connected to said recording head and each of the remaining signal components being produced at one of said playback heads;

channel means for conveying each time varying Fourier frequency component to an inverse Fourier transformer, said transformer comprising second analog means for combining the continuously time varying Fourier components C,,,(t) together so as to reconstitute the original time varying function fit);

said combining means consisting of a network of operational amplifiers. it 22. Apparatus as set forth in claim 21 wherein said inverse Fourier transformer consists of a network of operational amplifiers and attenuators. 

1. In a method for continuous analog analysis for a time varying input function f(t) in the time domain into at least one of its Fourier Frequency components Cm(t) the steps of: continuously producing N signal components fm(t)/N ranging from f0(t)/N to fN 1(t)/N wherein each signal component except f0(t)/N is delayed an amount T from the previous signal component; continuously combining separated ones of said signal componeNts (fm(t)) fm(t)/N together so as to produce at least one continuous output function Cm(t) corresponding to one of the Fourier frequency components of the time varying input function f(t) to be frequency analyzed.
 2. The method of claim 1 wherein the step of combining said separated components fm(t)/N together so as to produce at least one output function Cm(t) corresponding to one of the Fourier frequency components of the time varying input signal f(t) comprises the step of, combining the time delayed continuous signal components fm(t)/N together according to at least one of the component products of the matrix product: DF(t) where D is the discrete Fourier transform matrix and F(t) is the delayed input signal component matrix.
 3. The method of claim 2 wherein the step of combining the delayed signal components fm(t)/N according to at least one component product of the matrix product includes the step of, adding the delayed signal components fm(t)/N according to at least one of the component products of said matrix product.
 4. The method of claim 2 wherein D is defined by the matrix:
 5. The method of claim 1 wherein N is a power of
 2. 6. The method for reconstituting a time varying function f(t) in the time domain from a plurality of N number of its continuously time varying Cm(t) Fourier frequency components in the frequency domain the steps of: continuously combining the separate Fourier frequency components Cm(t) according to the product of the matrix equation f(t) DC(t) where D is the discrete Fourier transform matrix and C(t) is the matrix for the continuous separate Fourier frequency components.
 7. The method of claim 6 wherein the step of combining a plurality of the continuous time varying Fourier frequency components Cm(t) together so as to reconstitute the original time varying function f(t) includes, combining the separate Fourier frequency components Cm(t) according to the product of the matrix Equation: f(t) DC(t) where D is the discrete Fourier transform matrix and C(t) is the matrix for the continuous separate Fourier frequency components.
 8. The method of claim 7 wherein the D matrix is defined by the matrix equation:
 9. The method of claim 8 where
 10. In an analog apparatus for continuous frequency analysis of a time varying input function f(t) in the time domain into at least one of its Fourier frequency components Cm(t): means for continuously generating N signal components fm(t)/N ranging from f0(t)/N to fN 1(t)/N wherein each signaL component except f0(t)/N is delayed an amount T from the previous signal component; means for continuously combining said signal components fm(t)/N together so as to produce at least one continuous output function Cm(t) corresponding to one of the Fourier frequency components of the time varying input function f(t) to be frequency analyzed.
 11. The analog apparatus of claim 10 wherein the means for combining said separated signal components fm(t)/N together so as to produce at least one output function Cm(t) corresponding to one of the Fourier frequency components of the time varying input signal f(t) includes, means for combining the time delayed continuous signal components fm(t)/N according to at least one of the component products of the matrix product: DF(t) where D is the discrete Fourier transform matrix and F(t) is the separated input signal component matrix.
 12. The analog apparatus of claim 11 wherein said means for combining the delayed signal components according to at least one component product of the matrix product includes, means for adding the delayed signal components according to at least one of the component products of said matrix product.
 13. The analog apparatus of claim 11 where D is defined by the matrix:
 14. The analog apparatus of claim 10 wherein N is a power of
 2. 15. Apparatus comprising: first analog means for continuously producing N time varying Fourier frequency components Cm(t) in the frequency domain of a time varying function f(t) in the time domain; said first analog means comprising: means for continuously generating N signal components fm(t)/N ranging from f0(t)/N to fN 1(t)/N wherein each signal component except f0(t)/N is delayed an amount T from the previous signal component; and means for continuously combining said signal components fm(t)/N together so as to produce N varying Fourier frequency components Cm(t); channel means for conveying each time varying Fourier frequency component to an inverse Fourier transformer, said transformer comprising second analog means for combining the continuously time varying Fourier components Cm(t) together so as to reconstitute the original time varying function f(t).
 16. The apparatus of claim 15 wherein said means for combining a plurality of the continuously time varying frequency components Cm(t) together so as to reconstitute the original time varying function f(t) includes, means for combining the separate frequency components Cm(t) according to the component products of the matrix product DC(t) where D is the discrete Fourier transform matrix and C(t) is the matrix for the continuous separate Fourier frequency components.
 17. The apparatus of claim 16 wherein:
 18. The apparatus of claim 17 wherEin:
 19. In an analog apparatus for continuous frequency analysis of a time varying input function f(t) in the time domain into at least one of its Fourier frequency components Cm(t): means for continuously generating N signal components fm(t)/N ranging from f0(t)/N to fN 1(t)/N wherein each signal component except f0(t)/N is delayed an amount T from the previous signal component; said generating means including a tape recorder with a recording head and (N-2) staggered playback heads, said signal component f0(t)/N being connected to said recording head and each of the remaining signal components being produced at one of said playback heads; means for continuously combining said signal components fm(t)/N together so as to produce at least one continuous output function Cm(t) corresponding to one of the Fourier frequency components of the time varying input function f(t) to be frequency analyzed.
 20. An analog apparatus as set forth in claim 19 wherein said means for continuously combining signal components consists of a network of operational amplifiers.
 21. Apparatus comprising: first analog means for continuously producing N time varying Fourier frequency components Cm(t) in the frequency domain of a time varying function f(t) in the time domain; said first analog means comprising: means for continuously generating N signal components fm(t)/N ranging from f0(t)/N to fN 1(t)/N wherein each signal except f0(t)/N is delayed an amount T from the previous signal component; and means for continuously combining said signal components fm(t)/N together so as to produce N varying Fourier frequency components Cm(t); Said generating means including a tape recorder with a recording head and (N-2) staggered playback heads, said signal component f0(t)/N being connected to said recording head and each of the remaining signal components being produced at one of said playback heads; channel means for conveying each time varying Fourier frequency component to an inverse Fourier transformer, said transformer comprising second analog means for combining the continuously time varying Fourier components Cm(t) together so as to reconstitute the original time varying function f(t); said combining means consisting of a network of operational amplifiers.
 22. Apparatus as set forth in claim 21 wherein said inverse Fourier transformer consists of a network of operational amplifiers and attenuators. 